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Abstract 

We show that QCD Minkowski observables such as the e + e~ i?-ratio and the 
hadronic tau decay R T are completely determined by the effective charge (EC) beta- 
function, p(x), corresponding to the Euclidean QCD vacuum polarization Adler D- 
function, together with the next-to-leading order (NLO) perturbative coefficient of 
D. An efficient numerical algorithm is given for evaluating R, R T from a weighted 
contour integration of D(se lS ) around a circle in the complex squared energy s-plane, 
with p{x) used to evolve in s around the contour. The EC beta-function can be 
truncated at next-to-NLO (NNLO) using the known exact perturbative calculation 
or the uncalculated N 3 LO and higher terms can be approximated by the portion 
containing the highest power of b, the first QCD beta-function coefficient. The dif- 
ference between the R, R T constructed using the NNLO and "leading-b" resummed 
versions of p(x) provides an estimate of the uncertainty due to the uncalculated 
higher order corrections. Simple numerical parametrizations are given to facilitate 
these fits. For R T we estimate an uncertainty 5a s {m^) ~ 0.01, corresponding to 
<5a s (M|) ~ 0.002. This encouragingly small uncertainty is much less than rather 
pessimistic estimates by other authors based on analogous all-orders resummations, 
which we demonstrate to be extremely dependent on the chosen renormalization 
scheme, and hence misleading. 



PACS: 12.38.-t; 12.38.Bx; 12.38.Cy; 13.35.Dx. 
Keywords: QCD, resummation, tau decay, a s . 



1 Introduction 



There has been extensive recent interest [1-3] in the possibility of using measurements 
of R T , the total hadronic decay width of the r lepton normalized to the leptonic decay 
width, for precise determination of the renormalized strong coupling a s (M§) (or more 
fundamentally Aj^g). For this purpose R T apparently possesses a number of advantages 
compared with other QCD observables. It is an inclusive quantity which can be computed 
in QCD using the operator product expansion (OPE) supplemented by analyticity [4- 
7]. It has been calculated in QCD perturbation theory to next-to-next-to-leading order 
(NNLO) O(af) H 0- Power corrections are expected to be small || §], and since the r 
mass is below the threshold for charmed hadron production only the light quarks u, d, 
s are active, so QCD with Nf = 3 massless quark flavours should be applicable. R T can 
be rather accurately determined from the measured electronic branching ratio of the r or 
from the r lifetime |TIJ. In evolving up in energy scale from a s (m^) to a s (M§), which 
is customarily quoted in global comparisons, the fractional error in a s (Aj^g) is reduced. 
Measurement of the hadronic width of the Z° to directly determine a s (M|) shares the 
same advantages of being inclusive, calculated to NNLO in perturbation theory and having 
small power corrections, but suffers significant corrections from heavy quark masses, and 
much larger systematic experimental errors. 

Despite these undoubted advantages possessed by R T as a means of determining a s , the 
relatively low energy scale involved, s = m%, might lead one to expect sizeable corrections 
from uncalculated 0(ocf) and higher orders in perturbation theory. To assess the effect 
of these terms with our present limited state of knowledge one can employ a, necessarily 
approximate, all-orders resummation of the QCD perturbation series. A well-motivated 



framework for this is provided by the leading-6 approximation [11, 12], also sometimes 
referred to as naive non-abelianization |13|, |I4 |. In this approach the portion of perturba- 
tive coeffcients containing the highest power of b = g(lliV-2iV/), the first beta-function 
coefficient for SU(A^) QCD with Nf active quark flavours, is resummmed to all-orders. 
This can be accomplished given exact large- Nf all-orders results [15-17]. This technique 
has been applied to the QCD vacuum polarization Adler D-function [|i~5| , |i~6f , and its 
Minkowski continuations, the e + e~ i?-ratio and R T . Deep inelastic scattering (DIS) sum 
rules [nj] and heavy quark decay widths and pole masses have also been discussed [[L8|, |i~9|[ . 

In several recent papers it has been claimed that applying the leading-6 resummation 
to R T indicates rather large perturbative uncertainties |L8], ^(J, [H]. Indeed the estimated 
uncertainty in a s (M|) is of the same order as that normally quoted in determinations 
from jet observables at LEP and SLD [p2| . 

In a recent paper [23] we have pointed out that a straightforward resummation of the 



leading-6 terms of the kind employed in references |L8], pl| , is renormalization scheme 
(RS) dependent. This occurs because the compensation mechanism between the renor- 
malization group (RG) improved coupling and the perturbative coefficients is destroyed 
by retaining only the leading-6 terms. As a result the 'naive' leading-6 resummation is not 
RS-invariant under the full QCD renormalization group (RG). Whilst at large energies 
the resulting ambiguities are mild, at s = this RS dependence is serious and in our 
view invalidates the rather pessimistic conclusions of these references regarding the likely 
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uncertainty in a s (M^) determined from R T . 

In reference |23f we proposed an improved RS-invariant resummation based on ap- 
proximating the unknown effective charge (EC) beta-function coefficients by the portion 
containing the highest power of b. Approximated perturbative coefficients in any RS can 
then be obtained using the exact QCD RG. The leading-6 effective charge beta-function 
can be resummed using exact all-orders large- Nf results. 

The difference between the exact NNLO result for R T in the effective charge RS, 
and the RS-invariant all-orders resummation indicates a rather small uncertainty due 
to the approximated higher order terms, and the estimated uncertainty in a s (M§) is 
correspondingly small, 5a s (M§) ~ 0.001. 

In this paper we wish to explore the perturbative uncertainty in R T in somewhat more 
detail. Both R T and e + e~ .R-ratio can be represented by a contour integral involving 
D(se l9 ), where D(—s) is the Euclidean Adler .D-function, around a circle, cut along the 
positive real axis, 9 = —tt to 9 = n, in the complex s-plane Q]. Here s = m 2 for R T . 
Conventional perturbation theory involves an expansion in a s (s) obtained by re-expressing 
a s (se ld ) as an expansion in a s (s) which is then truncated. Alternatively one can simply 
numerically perform the contour integration for the a k s {se ie ) terms up to a given order 
[T[]. This procedure includes in addition to conventional fixed-order perturbation theory 
a resummation of analytical continuation terms. A subset of these terms involve powers 
of the first beta- function coefficient, b, together with ir 2 factors, and are resummmed to 
all-orders in the leading-6 approach. In addition, however, there are potentially large 
contributions involving higher beta-function coefficients |24|. It would seem sensible, 



therefore to perform the improved RS-invariant resummation for D(se ze ), and numerically 
evaluate the contour integral. In this way additional analytical continuation terms not 
captured in the leading-6 resummation are included exactly. This can then be compared 
with the exact NNLO result for D(se 10 ) in the EC scheme, with the contour integral 
again numerically evaluated. Since in both cases the analytical continuation terms are 
resummed the difference should be indicative of the effect of the approximated higher 
effective charge RS invariants for D beyond NNLO. 

The plan of the paper is as follows. In section 2 we shall intoduce the contour integral 
representations of the e + e~ i?-ratio and R T in terms of the Adler .D-function. Using 
Taylor's theorem we can then expand R and R T in terms of D(s) and its energy derivatives, 
which in turn can be expressed in terms of the effective charge beta-function for D and 
its derivatives. These results can be easily used to express the perturbative coefficients of 
the Minkowski quantities R and R T in terms of those of the Euclidean Adler .D-function 
and its effective charge RS invariants. We have compared these with existing expressions 



available in the literature [24]. We also derive relations between the EC invariants for 
R and R T and these for D. In section 3 we briefly review the basis of the RS-invariant 
resummation proposal |[23|| . The contour integrals for the R and R T are evaluated by 
using Taylor's theorem succesively to evaluate D at a series of values of complex s around 
the unit circle contour of integration. A Simpson's rule numerical integration is then 
performed. The translation of D in complex s involves the effective charge beta-function 
and its derivatives. This function can be truncated or its leading-6 terms resummed [p3 |. 
In section 4 fits to the experimental data for R and R T are performed to determine a s from 
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fixed-order and resummed perturbation theory. In section 5 we conclude by comparing 
the resulting values and estimates of the perturbative uncertainty with those suggested 
by other approaches. 



2 Contour integral representation of Minkowski ob- 
servables 

The two quantities with which we shall be concerned are defined as follows. 
The e + e~ .R-ratio is the observable 

_ <r(e + e~ — > hadrons) 
a(e+e- -+ ' 

In SU(iV) QCD perturbation theory 

R(s) = N £ Q) (l + -C F R) + Q f \ ~R , (2) 

where Q f denotes the electric charge of the quarks and the summation is over the flavours 
accessible at a given energy, s is the physical timelike Minkowski squared momentum 
transfer. The SU(iV) Casimirs are defined as C A = N,C F = (N 2 - 1)/2N. 

R denotes the perturbative corrections to the parton model result and has the formal 
expansion 

R(s) = a + na 2 + r 2 a 3 + ■■■ + r k a k+1 + • • • , (3) 

where a = a s (^ 2 )/ir denotes the renormalization group (RG) improved coupling. The MS 

scheme with /i 2 = s is often chosen. The R contribution first enters at 0(a 3 ) due to the 
existence of diagrams of the "light-by-light" type. 

The ratio R T is defined analogously using the total r hadronic width as 

„ r(r — ► u T + hadrons) , , 

r(r -> v T e v e ) 

Its perturbative expansion has the form 

5 a{m 2 ) 



R T = N(\V ud \ 2 + \V us \ 2 )S_ 



EW 



12 vr 



+ R T 



(5) 



where the V u d and V us are CKM mixing matrix elements, with |\4rf| 2 + |Kis| 2 ~ 1- Since 
the energy scale, s = m 2 , of the process lies below the threshold for charmed hadron 
production only three flavours u, d, s, are active. a{m 2 ) denotes the electromagnetic 
coupling p5j and Sew — 1-0194 p6| denotes further electroweak corrections. R T has a 



perturbative expansion of the form of equation (3) with coefficients which we shall denote 
rj.. There are no "light-by-light" corrections for R T since (SQ/) 2 = for u, d, s active 
quark flavours. 
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These two Minkowski quantities can both be expressed in terms of the transverse part 
of the correlator H(s) of two vector currents in the Euclidean region, 



n(s)(<7 M ?„ - g,uq 2 ) = 4vr 2 * / rf 4 xe^(0|T{^(a:)>(0)}|0) 



(6) 



where s = —q 2 > 0. In order to avoid an unspecified constant, it is convenient to consider 
the related Adler Z)-function, 

D(s) = -s^-n(s) . (7) 

In perturbation theory D can be written in the form of equation (2) involving perturbative 
corrections D with an expansion as in equation (3) involving coefficients dk, and "light- 
by-light" corrections D. 

R, R T and related Minkowski quantities such as spectral moments |J can all be written 
in terms of a weighted contour integral of D{se ld ) around a unit circle in the complex s- 
plane [g. 

Denoting such a generic Minkowski observable as R we have 



R(so) = tt d6W(6)D(s e l ") 
In 



(8) 



where W(8) is a weight function which depends on the observable R. For R we have 
W{6) = 1, and for R T 

W T (9) = (1 + 2e ie - 2e 3ie - e 4ld ) , (9) 

and So = m 2 .. 

A novel representation for R can be obtained by using Taylor's theorem to expand 
D(se td ) around s = sq. This yields 



1 f7r ( oo 

R(s ) = - J^d9W(9) \ D(s ) + 



,0' 



n=l 



n\ d In s n 



D(s) 



(10) 



The derivatives in equation (10) can be recast in terms of the effective charge beta-function 
p(D) |27], |28| , and its derivatives. p(D) is defined by 



dD{s) 
dins 



6,- 



(D 2 + cD 6 + p 2 ,D 4 + • • • + p k D 



k+2 



(11) 



Here b = ^(HN-2Nf) is the first coefficient of the beta-function, and 



7C| 11 C A C F 5^ 



8 b 



+ l c A + l c F 



(12) 



is the second universal beta-function coefficient. The higher coefficients p2, P3, ■ ■ -, in equa- 
tion (11) are RS-invariants and may be expressed in terms of the perturbative coefficients 
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of D, dk, together with the beta- function coefficients, cj., which define the renormalization 
scheme employed to define the RG improved coupling a ||29|| . Thus 



da(p 2 ) 
d In p 2 



-- (a 2 + ca 3 + c 2 a 4 + 



+ c k a k+2 + •••)• 



(13) 



The effective charge (EC) scheme corresponds to the choice of coupling D = a. The first 
two EC invariants are given by 



Pi = c 2 + d 2 — cdi - d\ 

p 3 = c 3 + 2d 3 - Ad ± d 2 - 2d x p 2 - cd\ + 2d\ . 



(14) 



We note that in references |23], |28|, to which the reader is referred for additional discussion 
of the EC beta-function, the dependent energy variable was taken to be Q, whereas we 
are employing s = Q 2 in this paper, hence there are additional factors of ~ in equations 
(11) and (13). 

Using equation (11) one can then show that the energy derivatives in equation (10) 
can be rewritten as 



dlns n 



D(s) 



p(x) 



dx 



1 n-l 

p{x) 



n > 



(15) 



x=D(s ) 



Thus finally we can write equation (10) in the form 



-%b \ w, 



71 

n\ 



p{x) 



R(s ) = D(s ) + J2 

Here w n denotes moments of the weight function W(8), 

d9 n W{6) 



dx 



-I n-l 

p(x) 



(16) 



x=D(sq) 



1 

2^ 



(17) 



For R setting W{6) = 1 yields w n = ix n /{n + 1) (n even), w n = {n odd). The first two 
terms in the sum of equation (16) are then 



R(s ) = D(s ) 



7l 2 b 2 , 7T 4 6 4 ,, 

-pp' + 



24 



1920 



p^ + App'p" 



P 2 P"') + ■ ■ 



Primes denote differentiation of p(x) with respect to x, evaluated at x = D(s ). Succes- 
sive terms are RS-invariants resulting from the resummation to all-orders of analytical 
continuation terms proportional to ir 2 b 2 , 7r 4 6 4 , • • •, respectively. 



For R T the weight function W T (8) of equation (9) has moments w\ 
. Then equation (16) yields 



19? 



12 ' W 2 



265 
3 72 ' 



R T = D{m 



2N 196 /tt 2 265 U 2 . 
-I n — I I — pp H 



24 



P 



72 J 8 



(19) 



/,From equation (16) we see that Minkowski observables are naturally expressed in terms 
of the Euclidean Adler Z)-function and its EC beta-function. Given an all-orders definition 
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of p{x) one can discuss the radius of convergence of the sum in equation (16) in D(s ). 
This is an interesting question which could be directly addressed using the leading-6 



resummmation of p on which the RS-invariant resummations of reference |23[ are based. 
However, one can anticipate a rather restricted radius of convergence by making the one- 
loop approximation p(x) = x 2 . One then finds 

~ 2 MD. , , 

-Ronoioop = — arctan(— — ) , (20) 
on 2 

for the result of resumming the analytical continuation terms which only involve the 
lowest beta-function coefficient. This suggests that the radius of convergence is limited 
by D < ^- [|I], |30|]. For Nf = 3 this gives a radius of convergence 5- ~ 0.141 • • •, which is to 
be compared with D(m 2 ) ~ 0.1 So that the expansion will not be useful for evaluating R T 
using the leading-6 resummation of p(x). As we shall discuss in the next section, however, 
we shall be able to make use of the Taylor's theorem approach to evaluate D(se ie ) at 
closely spaced intervals around the integration contour using the resummed p(x). 

To conclude this section we note that the expansion of equation (16) is of use in 
straightforwardly relating the R, R T Minkowski perturbative coefficients rj., r k to the 
Euclidean c4 coefficients of D. One simply expands equation (16) as a power series in D 
then on substituting D = a + d\a 2 + d 2 a 3 + • ■ •, and isolating the coefficient of a k+1 one can 
directly obtain r k , r k in terms of di<k, Pi<k and c. The resulting calculated expressions are 



in agreement with the results available in the literature for k < 5 ||24j| , on using equations 
(14) to re-express beta-function coefficients c% in terms of pi invariants. 

In clarifying the connection between the various versions of fixed-order perturbation 
theory to be compared in section 4 it will be useful to relate the EC RS-invariants p k 
and p k T corresponding to R and R T to the p k (hitherto pp.) connected with D f31|| . This 



can easily be accomplished by first evaluating r k and r T k in the EC scheme for D, so that 
d\ = d 2 = ■ ■ ■ = 0. These r^^di = 0) and r T k {di = 0) are simply the coefficient of D k+1 on 
the right-hand side of equation (16). One can then use the analogue of equation (14) for 
R, R T with Ck = p k to obtain the required relations. Expressions for p k and p k T (k < 6) 
are included in Appendix A. 



3 Numerical evaluation of the contour integral 

In this section we shall reformulate the Taylor's theorem expansion approach of the last 
section to obtain a tractable method for numerically evaluating the contour integral, ap- 
propriate not only when D is truncated at fixed-order in perturbation theory but crucially 
also allowing us to perform the RS-invariant all-orders resummation Z)( L *) of reference |p3[| . 

For ease of reference we shall begin by briefly reviewing the leading-6 resummations. 
The reader is referred to reference [^] for full details. 



For a wide range of "quark-initiated" [14] QCD observables the perturbative coeffi- 
cients can be organized as polynomials in the number of quark flavours, Nf. That is 
assuming such an observable D(s) with the perturbation series 

D(s) = a + dia 2 + d 2 a 3 H h d k a k+1 H , (21) 
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we can write 

4 = efWj? + d?- 11 ^ 1 + • • ■ + 4° ] • (22) 
By substituting Nf=(Y-N — 36) equation (22) can be recast in powers of b. 

4 = d[ k) b k + df-% k - 1 + • ■ ■ + 4 0) • (23) 

Since 4^ = ( — l/3) fe 4 ' exact all-orders large- Nf results can be used to perform a 
"leading-6" resummation, 

oo 

= £ d ( k L) a k+1 , (24) 

k=0 

where d[ L) = d[ k) b k (d L) = 1). 

may be defined as the principal value (PV) regulated Borel sum 

/•oo 

Z)( L )(a)=PV/ dze- z/a B[D^](z) , (25) 
Jo 

B[D^](z) denotes the Borel transform, which potentially involves an infinite set of poles 
at z — Zi = j (I = 1, 2, 3, • • •) corresponding to infra-red renormalons (IRz), and at z=— z\ 
corresponding to ultra-violet renormalons (UV;). In the specific case of the Adler D- 
function, D^ L \ IR4 is not present reflecting the absence of a relevant operator of dimension 
two in the operator product expansion (OPE) for vacuum polarization [jTT], |l6fl . IR 2 is a 
single pole and the remaining singularities are double poles. Expressions for the residues 
are given in reference [jT^]. It is then straightforward to evaluate equation (25) in terms 
of the exponential integral functions, 

/OO f>~ t 

dt— , (26) 
-X t 

where for x > the Cauchy principal value is taken. The UV; singularities may then be 
expressed in terms of Ei(— Fzi), where F = -, and the IR/ singularities involve Ei(Fzi). 
Corresponding expressions for D^ L \F)\uy and D^ L \F)\m are given in equations (48) and 
(49) respectively of reference |T2]| . 

As pointed out in reference p3| the resummation of equation (24) is ambigu- 
ous due to its RS-dependence. In particular if, as in the case of D, the exact NLO and 
NNLO coefficients are available it would seem sensible to include them and approximate 
only the unknown 4? d^, ■ ■ •, higher coefficients by d^\ Unfortunately, however, the 
resummed result explicitly depends on the RS chosen for evaluating the exact di,d 2 co- 
efficients. This is analogous to the ambiguity encountered in matching leading logarithm 



resummations of jet observables to exact fixed-order perturbative results p2 |. In both 



cases the difficulty may be avoided by performing a resummation of the EC beta-function 



23 , f28f . The unknown N 3 LO and higher EC beta-function coefficients p 3 , p 4 , • ■ -, in equa- 



tion (11) are approximated by retaining only the portion involving the highest power of 
b, p^ = p^b k . The p^ can be obtained to all-orders from the large- Nf result for 4^- 
If the NNLO invariant p 2 is known exactly it can be included. One then arrives at the 
resummed EC beta-function 

00 

p i - L *\x) =x 2 (l + cx + p 2 x 2 + Pk ] x k ) • (27) 

k=3 
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The improved RS-invariant resummation D^ L *\s) can then be obtained as the solution 
of the integrated beta-function equation 



£)0) 



-t-cln 



--In 
1 + cD( L *) 2 A 



MS 



D (L*) 

dx 



p( L *)(x) 



x 2 (l + cx) 



The exact NLO coefficient d\ occurs in the RS-invariant combination 

b 



Po 





2 Ars 



-d?V = a ). 



(2? 



(29) 



The convention assumed for A^ in equation (28) differs from the standard definition 
by the iVj-dependent factor (2c/o)~ c / b |p9|| . If D^ L *\s) is expanded in the coupling a 
appropriate to some RS one then obtains 



a 



{dr = i) 



(30) 



k=0 



where now d[ L *^ = d\ and d 2 



(L*) 



d 2 reproduce the known coefficients and the approx- 
imated o^*'' and higher coefficients may be obtained in any RS from the approximated 
invariants using the exact QCD RG. For instance if we label the RS by d x and the 
beta-function coefficients c 2 , c 3 , • • •, appearing in equation (13) we have [|3[ 

dt\di, c 2 , c 3 ) = d\ + b -cd\ + (3p 2 - 2c 2 )d x + i(p< L) - c 3 ) . 



(31) 



This differs from the exact d 3 only in the unknown p 3 term (we note in passing that c 3 
has now been computed the known p 2 has been exactly included. In this approach 
the maximum available exact information is included in an RS-invariant manner. 

It finally remains to perform the resummation p^ L *^ using the explicit expressions for 
D^(F). Defining 



p W(*)^ 2 (i + £p1 z V) 

k=2 



(32) 



and using the chain rule to relate the beta-function in two different RS's one has 

x 2 dD( L \a 



p( L )( x ) = [a (L) (x) 



do 



(33) 



a=a( L ~) (x) 



where aS L \x) is the inverse function to D^(a), i.e. D^ L '(a^(x)) = x. 



p( L \x) 



can then be straightforwardly defined from D^ L \F). The first step is to nu- 



merically solve 

D {L \F(x)) = x , 
to obtain F(x). Recalling that F = - one can then determine 



(34) 



(35) 



F=F(x) 
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by differentiating the explicit D^ L \F) expressions in terms of Ei functions. Finally on 
comparing equations (32) and (27) one has 

( x ) = pM (x) + cx 3 + p [ 2 NL) x 4 , (36) 

where = p 2 — p% . p^ L *\x) can then be inserted in equation (28) and the inte- 

gration performed numerically. On solving the transcendental equation the RS-invariant 
resummation can be evaluated. 

We now turn to the problem of evaluating the improved resummations ES L **^ and 
jl(L**) j22j w h ere contour integration of equation (8) is performed with D^ L *\se %e ) in 
the integrand, 

R( L **\s ) = d9W(9)D {L *\s e ie ) . (37) 

27T J -it 

To perform the contour integration numerically one can split the range from 9 = 0, ir into 
K steps of size A9 = j^, and perform a sum over the integrand evaluated at 9 n = nA9 
n — 0, 1, • • • , K. So that 

R(s ) ~ — \W(0)D(s ) + 2Re ^ W(9 n )D(s n )] , (38) 
27r L „=i J 

where s n = soe inAd . In practice we perform a Simpson's Rule evaluation. 

An efficient strategy |33J] is to start with the exact D(s ) and evolve D(s n ) to D(s n+ i) 
using Taylor's theorem. Thus defining x n = D(s n ) we have 

Xn+1 = Xn . i^ bpM _ 

2 o 

+ ! ^^& 3 |p(x„)(p / (x n )) 2 + (P(X„))VW| + O(A^) 4 + • • • , (39) 

analogous to equation (16). If equation (39) is truncated by retaining its first m terms 
one anticipates an error ~ K ^- 2 in equation (38). 

To evaluate R T to four significant figure accuracy retaining the first four terms in 
equation (39) we required 100 steps. 

The method obviously can also be used to evaluate the contour integral when D is 
represented by fixed-order perturbation theory in the coupling a(soe ), 

D(s e ie ) = a(s e ie ) + rfia 2 (s e ie ) + d 2 a 3 (s e ie ) . (40) 

One can start with a(so) and evolve a(s n ) to a(s n+ i) using equation (39) with p(x) replaced 
by the truncated beta-function in the corresponding RS, 

B{x) = x 2 + cx 3 + c 2 x 4 . (41) 

In standard approaches [TT| the contour integral is performed by solving the integrated 
beta-function equation with complex renormalization scale s n for a(s n ) at each integration 
step, and takes much longer to evaluate. Reference |HJ considers in some detail the RS 
dependence of the contour integral. 
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Determining PS L **"> with this approach is now relatively straightforward. For some 
given Ajyjg one evaluates L)( L *\so) [^3] as we have reviewed. The truncated equation 
(39) is then used to obtain x± = D^ L *\si). This requires p( L *\x ) and some number of 
derivatives. p( L >(xo) can be determined given the D^ L \F) expressions of reference [[Cql , 
and using the numerical inversion route of equations (34) and (35). p^', p^ L '", • • • can 
then be obtained by successive differentiation of equation (34) with respect to x. One 
finds 



and 



pV»(x) = ^ (L)//(F(X)) ) 2 - ly ' " V ^ y) (43) 



where primes denote differentiation with respect to F. Thus, once F(x) has been deter- 
mined from equation (34) no further transcendental equations need to be solved and the 
explicit expressions for D^ L \F) can be repeatedly differentiated to obtain p^', p^ L '", • • •■ 
Finally p^ L *\ p^ L *^', p^*- 1 ", ■ • •, can be obtained using equation (36) and its derivatives. 
For instance 

p^\x) = p (L) '(x) + 3cx 2 + 4p ( 2 NL) x 3 . (44) 

The only remaining difficulty is that x\ is now complex, and so at subsequent steps it is 
unclear how to obtain p( L *'(x n ), since D^ L \F) is only defined for real F. One needs to 
replace the Ei(x) defined in equation (26) by the generalized exponential integral functions 
Ei(n, w) for complex w, used to evaluate and RljP' in reference fl2[]. These are defined 
for Rew > by 

/oo e~ wt 
dt ~^r • ( 45 ) 

For Rew < they are defined by analytical continuation to arrive at a function analytic 
everywhere in the cut complex w-plane except at w — 0, and with a branch cut running 
along the negative real axis. 

To define D^ L \F) correctly for complex F one needs to replace Ei(— Fzi) in equation 
(48) of reference JL2] for D (L \F)\ VV by -Ei(l, Fzj). In equation (49) of reference Jl2 



for D^ L \F)\ m one replaces Ei(F^) by — Ei(l, —Fz{) + ivrsign(ImF^). In this way as F 
becomes real one avoids ±z7r contributions from the discontinuity across the branch cut 
along the negative real axis and re-obtains D^ L '(F) for real argument. 

With D^ L \F) re-defined for complex arguments in this way xi,X2, • • •, can be suc- 
cessively obtained. At each step one needs to solve the complex-valued transcendental 
equation 

L)M{F n )=x n , (46) 

and F n is then used to construct p^ L *\x n ) and its derivatives using equations (35), (42) 
and (43). The required computing time is dominated by that required for the solution of 
equation (46), and is comparable to that needed for the conventional approach in fixed- 
order perturbation theory, where the complex-valued integrated beta-function equation is 
numerically solved at each step. 
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Figure 1(a): Comparison of two versions of fixed-order EC perturbation theory "con- 
tour-improved" M L **W(EC) ("+") and R^ L **^ n \EC) ("x") with the RS-invariant re- 
summation R( l **\sq) (dashed line) at ^s = 91 GeV. 



We have checked that evaluating the contour integral with and p^ reproduces 
values in numerical agreement with the R> L > and R^' expressions of reference |12| . 



In the next section we shall compare the "contour-improved" RS-invariant resumma- 
tions and R^ L **\ with "contour- improved" fixed-order results obtained by trun- 

cating £)( L *) at n th order in the EC scheme, that is in equation (39) p(x) is taken to 
be the truncation of p^ L *'{x) in equation (27), retaining terms up to x n , and the input 
(s ) is obtained by solving equation (28) with the truncated p^ L *\x). We shall 
denote these by M L **^{EC) and R[ L **^(EC) for n > 3, and for n = 1,2 where the 
exact p k are used by R^(EC) and R^{EC). 

These "contour-improved" evaluations are to be compared with conventional fixed- 
order perturbative truncations R^ L **^ n \EC) and R^**^ n \EC) obtained by integrating 
up the n* h -order truncated EC beta-functions p R > R r( L **)^ with coefficients p R ' Rr ( L **) 0D _ 
tained using equations (51), (52) in Appendix A, with the exact p® and using pj?^ for 
k > 2. By truncating p R ' Rr one omits an infinite set of exactly-known and numerically 
important analytical continuation terms which are included in the "contour-improved" 
resummations. 



4 Numerical results 



In Figures l(a)-(c), for y/s = 91,5,1.5 GeV, respectively, we compare the "contour- 
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Figure 1(b): As for Figure 1(a) except at y/s = 5 GeV. 
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Figure 1(c): As for Figure 1(a) except at y/s = 1.5 GeV 
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Figure 2(a): Comparison of fixed-order EC perturbation theory D^ L *^ n \EC) ("+") with 
the RS-invariant resummation L)^ l *\sq) (dashed line) at y/so = 91 GeV. 

improved" resummation R^ L **\so) with the two versions of fixed-order perturbation the- 
ory, "contour-improved" R( L **W(EC) and R^ L **^ n \EC) (for n > 2), described in the 
last section. Values of A§L = 200 MeV, = 279 MeV and A§L = 320 MeV are used. 
These assume flavour thresholds at m b = 4.5 GeV and m c = 1.25 GeV. 

As can be seen at all energies and in low orders the "contour-improved" fixed-order 
results (denoted "+") are significantly closer to the resummation i?( L **) (horizontal line) 
than the conventional fixed-order results (denoted "x"). This is completely understand- 
able since both the RS-invariant resummations and the contour-improved fixed-order re- 
sults sum to all-orders known analytical continuation terms, as discussed above. The 
unnecessary truncation of these terms evidently greatly worsens the performance of n = 2 
NNLO fixed-order perturbation theory, whilst in higher orders both versions of fixed-order 
perturbation theory approach each other, and both track the RS-invariant resummation. 
Eventually, of course, both versions will breakdown as the leading UVi renormalon singu- 
larity asserts itself. Since n = 2 represents the highest order for which exact calculations 
exist at present, "contour-improvement" is clearly essential if reliable NNLO determina- 
tions of a s [M\) are to be made. 

In Figures 2(a)-(c) we plot D^ L *\s ) (dashed line) and D (L *^ n \EC) (denoted "+") 
at yfs = 91,5,1.5 GeV, respectively. These represent the input values of D(sq) fed 
into the contour integration to produce the plots in Figures l(a)-(c). We note that the 
fixed-order results in Figures 2 show a clear oscillation above and below the resummed 
result. This is a reflection of the alternating sign factorial behaviour contributed by 
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Figure 2(b): As for Figure 2(a) except at y/s = 5 GeV. 



0.140 



.135 



.130 



. 125 



. 120 



.115 




Figure 2(c): As for Figure 2(a) except at y/s = 1.5 GeV. 
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Figure 3: As for Figure 1(a) except for R T at a/s = 1.777 GeV. 



the leading UVi renormalon, which in the case of D is a double pole, in the leading- 
h approximation. A similar oscillatory behaviour is also evident for the conventional 
fixed-order perturbative approximants for R in Figures l(a)-(c), but with much smaller 
amplitude. This is because for R the UVi singularity is softened to a single-pole, again in 
the leading-6 approximation. As a result one expects r n /d n ~ - asymptotically [11], pfl , 
and correspondingly p^l Pn — ■ Notice that the "contour- improved" fixed-order results 
which partially resum higher-order contributions do not exhibit the simple oscillatory 
behaviour. 

In Figure 3 we give the analgous plot to Figures 1 for R T , assuming A^ = 320 MeV 
as before, y/s = 1.777 GeV. If we compare with Figure 1(c) for R at the comparable 
energy y/s Q = 1.5 GeV, we see a deterioration in the behaviour of both versions of fixed- 
order perturbation theory. The change of weight function from W(8) = 1 to W T = 
(1 + 2e ie - 2e 3ie - e 4ie ) leads to much less convergent analytical continuation terms and 
the two versions of fixed-order perturbation theory no longer approach each other in higher 
orders. The contour-improved results are reasonably close to the resummation. Clearly 
"contour- improvement" is vital for reliable a s (m 2 T ) determinations. 

We now wish to use the difference between the "contour-improved" i?( L **) and Rw {EC) 
to estimate the uncertainty with which a s (M|) can be determined for the Minkowski ob- 
servables. Our main interest will be in R T which potentially gives the most accurate 
determination. To begin with, however, we consider R at ^/s = Mz (i.e. the hadronic de- 
cay width of the Z°). As in our fits in reference |3| we shall take R{M 2 Z ) = 0.040 ±0.004. 
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Figure 4: As for Figure 3 except R^**) is fitted to R e T xpt = 0.205 at y/s = 1.777 GeV. 



The fits to the three-loop NNLO MS a s (M 2 z ) are then a s {M 2 z ) = 0.122 ±0.012 from both 
R( L **1 and R®(EC). This is also the same result as obtained in reference 1231 using 
So at this high energy scale the "contour-improvement" has little effect. 

For R T we take _R^ ata = 3.64 ± 0.01 [ 2Tf . Correcting for the small estimated power 



corrections |20[ then yields R T = 0.205 ± 0.006. Fitting to the "contour-improved" RS- 
invariant resummation R^ L **^ then yields a s (m 2 ) = 0.339 ± 0.006, and fitting to the 
"contour-improved" NNLOEC result R^(EC) gives a s (m 2 T ) = 0.350 ± 0.008. Evolving 
the three-loop coupling a^ s using Bernreuther- Wetzel matching [|35| from Nf = 3 to 
N f = 5 with the flavour thresholds noted above yields a s (M z ) = 0.1214 ± 0.0007 and 
a s (M z ) = 0.1226 ± 0.0008, respectively. Using smaller quark masses (m c = 1.0 GeV 
and mi, = 4.1 GeV) at the bottom of the range quoted in |36[] to perform the evolution 
yields a s (M|) = 0.1222 ± 0.0007 and a s {M z ) = 0.1234 ± 0.0008, respectively. Choosing 
larger masses (m c = 1.6 GeV and nib = 4.5 GeV) at the top of the quoted range gives 
a s (M 2 ) = 0.1207 ± 0.0007 and a s (M 2 ) = 0.1218 ± 0.0008, respectively Thus one can 
estimate an uncertainty 6a s (M z ) ~ 0.002. 

In reference |23| we fitted to the same value of R T using the RS-invariant resummation 
which only includes analytical continuation terms at the leading-6 level, and found 
a s (m 2 ) = 0.328 ± 0.005, similarly fitting to NNLO EC fixed-order perturbation theory 
gave a s (m 2 ) = 0.320 ±0.005. In both cases the inclusion of exactly known analytical con- 
tinuation terms involving c, P2, ■ ■ ■ , via the "contour-improvement" serves to significantly 
increase the fitted a s (m^), and hence slightly increase a s {M z ). 
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In Figure 4 we repeat the plot of Figure 3 but with the increased value of = 429 
MeV, which results from fitting _R(. L **) to the data for R T . A marked deterioration in 
the performance of both versions of fixed-order perturbation theory is evident, although 
the "contour-improved" fixed-order results are still significantly closer to the resumma- 
tion. This serves as a warning that, at this low energy scale, relatively small changes in 
m 2 /(A^) 2 can produce significant changes in the accuracy of perturbation theory. 

In Figure 5 we extend the fits reported earlier. For 0.16 < R T < 0.25 we plot curves for 
the a s (m 2 ) obtained by fitting R^(EC) (dotted curve), R[ L **^ (solid curve), R^(EC) 
(dashed curve), to this value. 5a s {m 2 ) can then be estimated for given R T from the 
difference between the lower two "contour-improved" curves. Clearly 5a s (m 2 ) increases 
very rapidly as R T increases. We are very fortunate that apparently R T ~ 0.2, for which 
5a s {m 2 T ) ~ 0.01. 

We can compare the RS-invariant resummation R^ L **^ with the results obtained using 
a Le-Diberder Pich (LP) resummation [Q, that is evaluating the contour integral with 
D(m 2 T e i6 ) as in equation (40). Fitting R T = 0.205±0.006 to this £^(MS) yields a s {m 2 T ) = 
0.359±0.008 in reasonable accord with the NNLO EC "contour-improved" value a s (m 2 ) = 
0.350 ±0.008. 

Fitting to conventional NNLO fixed-order perturbation theory in the MS scheme 
-R[ 2 )(MS), with ix = m T , yields a s (m 2 ) = 0.342 ± 0.006. No special significance should 
be ascribed to the numerical coincidence that this is close to the RS-invariant resumma- 
tion fit, a s {m 2 ) = 0.339 ± 0.006. Crucially these MS fits are strongly dependent on the 
assumed RS. In Figure 6 we show the R T versus a s {m 2 ) plots with three choices of scale 
/i = 2m T , m T , \m T (labelled 1, 2, 3, respectively). We also plot the curves for fitting to a 
LP resummation based on the NNLO MS expansion of D(m 2 e ie ) in a(4m 2 e t6 ), a(m 2 e ie ) 
and a(\m 2 e %e ) (labelled 4, 5, 6, respectively). As can be seen the i?^(MS) curves for 
different scales are very widely separated. The scale dependence of the LP resummations 
is seen to be much reduced compared to conventional fixed-order MS perturbation theory, 
but is still significant. 

For convenience we now present simple numerical parametrizations for the contour- 
improved resummations R[ L **\ w 2 \EC) and LP (i.e. i?^(MS) based on an expansion 
in a{m 2 e %e )), in terms of a s {m 2 ). We stress that a s {m 2 ) denotes the 3-loop NNLO MS 
coupling with scale /i = m T . 

Given x = R T (data) the fitted a s {m 2 ) is parametrized by 

a s {m 2 T ) = ttx + A 2 x 2 + A 3 x 3 + A 4 x 4 . (47) 

The numerical coefficients Ai for the different "contour-improved" versions of perturbation 
theory are tabulated in Table 1. These coefficients give a s (m 2 ) to a numerical accuracy 
of three significant figures over the range 0.16 < R T < 0.25 covered in Figure 5. 

We also present reverse fits. Given x = a s (m 2 ) the different approximations for R T 
are parametrized by 

R T = -x + A 2 x 2 + A 3 x 3 + A 4 x* . (48) 

7T 

The numerical coefficients Ai are again tabulated in Table 1. R T is accurate to three 
significant figures over the range 0.29 < a s {m 2 ) < 0.41. 
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Figure 5: R T versus a s (m 2 T ) for (EC) (labelled 1), (labelled 2) and R^(EC) 

(labelled 3). 
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Figure 6: R T versus a s (m 2 T ) for i?(. 2 )(MS) and i?P](MS), with three choices of scale 
/i = 2m T ,m T , \m T \ labelled (1, 2, 3) and (4, 5, 6) respectively. 
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Perturbative 
approximation 


At 


A 


A 2 


A 3 


A 4 


A 2 


A 3 


A 4 




-13.23 


+38.87 


-47.44 


+.6146 


+.2099 


+1.352 


RM(EC) 


-13.61 


+39.75 


-36.71 


+.09813 


+4.504 


-7.414 




-13.92 


+45.55 


-52.19 


+.4319 


+2.121 


-3.823 



Table 1: Numerical coefficients A, and A, to parametrize a s (m 2 ), given x = R T , and 
R T , given x = a s (m 2 ), respectively, for perturbative approximations R(p**\ Rfp(EC) and 
RW(MS). 



Finally in this section we wish to examine the performance of a straightforward leading- 
b resummation for R r . To emphasise the associated RS ambiguity we shall evaluate it for 
three MS scales p = Xm T , where A = |, 1, 2, as before. We then evaluate 

Rl L \F(a)) + r\ {NL) a 2 + r T 2 {NL) a 3 . (49) 

Here a denotes the MS coupling a(\ 2 m 2 ). R[ L \F) is given by the explicit expressions in 
equations (69) and (70) of reference |[2j. F(a) = ^ — 6(lnA + |). r^ NL ^ = r\ — rJ^ L \ The 
extra terms ensure that at NLO and NNLO the known exact rj", r\~ are included in the 
resummation. 

Fitting equation (49) to R T = 0.205 as before yields a s (\m 2 T ) = 0.475, a s (m 2 T ) = 0.306, 
and a s (4m^) = 0.233, for the three choices of RS. Evolving these all to m T using the three- 
loop MS beta-function, which is presumably appropriate since we are including the exact 
fixed-order results to NNLO, one obtains a s (m 2 ) = 0.297, 0.306, 0.322, respectively. Even 
if we restrict the beta-function to the one-loop leading-6 level, as advocated in references 
14], [18|, |3, we obtain a s (m 2 T ) = 0.323, 0.306, 0.304. Only if we leave out the NL correction 
terms in equation (49) and perform a pure leading-6 resummation do we uniquely obtain 
a Am 2 ) = 0.305 for all three RS's. 



5 Discussion and conclusions 

The essential point which motivates our approach is that the basic ingredient out of which 
the Minkowski observables R are built is the EC beta-function p(D(s)) defined in equation 

12 

(11). Using equation (7) one can see that this is proportional to ^j-p — -U(s), where n(s) is 
the fundamental correlator of two vector currents in the Euclidean region defined in equa- 
tion (6). If one specifies p(x), then given the NLO perturbative coefficient c?5 kls (/i 2 = s ) 
and assuming some value of A^jg, D(so) can be obtained unambigously on solving equation 
(28). There is no scale dependence since D(sq) only involves the RS-invariant combina- 
tion po( s o) m equation (29). Using equation (16) R(s ) is then also uniquely specified 
given p(x), where in practice the infinite sum is performed by numerically evaluating the 
contour integral, using p(x) to evolve D(se td ) around the circular contour of integration, 
as described in section 3. 
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Of course, the function p(x) is not known exactly. ^From NNLO calculations all that 
is known is the first three terms in its power series expansion, 

p(x) = x 2 + cx 3 + p 2 x^ H . (50) 

The uncertainty in predicting R(sq) is then to be estimated by making some approxima- 
tions for the unknown higher order terms indicated by the ellipsis in equation (50). We 
have chosen to approximate pk by for k > 3. These leading-6 contributions exactly 
reproduce pk in the large- Nf limit, and for p 2 are a good approximation in the large- N 
(or Nf ~ 0) limit |23j. Comparing the predictions for R T constructed from the NNLO 
p(x) in equation (50) and the leading-6 resummation indicates a moderate uncertainty 
5a s (m 2 T ) ~ 0.01 for R T ~ 0.2, which evolves up to 5a s (M|) ~ 0.002 and a central value 
a s (M|) = 0.122 in line with other a s measurements, which indicate a global average 
a s [M\) = 0.118 ±0.005 |22 . 



Our reassuringly small uncertainty 5a s (m^) ~ 0.01 is in stark contrast to other more 
pessimistic claims in the literature. Application of straightforward leading-6 resummations 
compared to exact NNLO fixed-order perturbation theory leads to a claim of 5a s (m 2 .) ~ 
0.05 in reference As we showed in section 4, however, there is a matching problem 



if one wishes to include the exactly known NLO and NNLO coefficients. As a result the 
6a s (m%.) estimate depends strongly on the renormalization scale chosen. This difficulty is 
avoided in our RS-invariant resummation approach, and originally motivated it. 

In reference [37J an overall uncertainty of 5a s {m 2 -) ~ 0.06 is claimed. These authors use 



an LP resummation together with an acceleration technique applied to the perturbation 
series to lessen the influence of the leading UVi renormalon. The resulting uncertainty 
is dominated by the choice of renormalization scale p. As we have pointed out above the 
only uncertainty in R T is due to our lack of knowledge of the uncalculated RS-invariants 
p 3 ,p 4 , Thus there is no scale dependence ambiguity. Since it is p(x) which is am- 
biguous one could attempt to improve the convergence of this series. The corresponding 
Borel transform has a UVi renormalon and one could try to use acceleration methods. 
Crucially, however, the resulting uncertainties would have to do with real ambiguities 
associated with the singularities of the Borel transform of p(x) in the Borel plane, and 
would not involve the unphysical and irrelevant renormalization scale p. The same crit- 
icism applies to reference |3£| which uses similar techniques to assess the perturbative 
ambiguity in R. 

We therefore conclude that there is no reason to suppose that R T suffers from serious 
ambiguities due to N 3 LO and higher terms which have yet to be exactly calculated. The 
techniques on which existing claims to this effect have been based are all severely RS- 
dependent, and their conclusions can be modified at will by making different ad hoc 
choices of renormalization scale. 
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A Relations between RS-invariants for R T and 
the Adler D-function 

Below we present the analytical continuation terms that link the RS-invariants for the 
two Minkowski observables R, R T to those of the Euclidean Adler D-function. 

For R we can relate the Minkowski invariants to the Euclidean invariants in the fol- 
lowing manner. 

pi = Pi-\f- 2 
Pi = Pi-Y2 CbV 

Pi = pf-^(8p? + 7c 2 )6V + ^V (51) 
Pi = Pi-^pi + 20pfc + 3cW + ^cfeV 

Pi = Pi-^(17P? + 28pfc+13(pf) 2 + 12pfc)6V 

+ — (99pf + 137c 2 )6 4 tt 4 - — L_ 
720 V 12 1 20160 

For R T we have more complicated relations. 
pf = p° + I 2 b 2 
pf = p *> + 5c/ 2 6 2 + I 3 b 3 
pf = pi + (8pi + 7c 2 )I 2 b 2 + 7cI 3 b 3 + hb 4 

pf = p f + (12pf + 20p^c + 3c 3 )/ 2 6 2 + (12pf + 16c 2 )/ 3 6 3 

+i(83/ 4 + 28/ 2 2 )c6 4 + I 5 b 5 (52) 
9 

pf = pi + (17p° + 28p°cn 2 + 13(p°) 2 + 12p°c 2 )I 2 b 2 

+-(39pf + 99p^c + 30c 3 )/ 3 6 3 
2 

+— (612p£ + 1081c 2 )/ 4 6 4 + 4( 234 P? + 277c 2 )/ 2 6 4 
36 18 

+^(93/ 5 + 60/ 3 / 2 )c6 5 + I 6 b 6 

o 

where for convenience we have assigned 
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1819 
3456 

246779 169 2 1 4 
165888 ~ 864^ + 360^ 

269203 1819 2 
55296 ~~ 3456^ 

392305009 973531 2 1859 
21233664 ~ 442368 * + 46080 
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